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^SJ ' We calculate the transition temperature in ultranarrow superconducting wires as a function of 

resistance, wire width and applied magnetic field. We compare the results of first-order perturbation 
(j^ ', theory and the non-perturbative resummation technique developed by Oreg and Finkel'stein. The 

' latter technique is found to be superior as it is valid even in the strong disorder limit. In both cases, 

the predicted additional suppression of the transition temperature due to the reduced dimensionality 
is strongly dependent upon the boundary conditions used. When we use the correct (zero-gradient) 
boundary conditions, we find that theory and experiment are consistent, although more experimen- 
tal data is required to verify this systematically. We calculate the magnetic field dependence of 
C ■ the transition temperature for different wire widths and resistances in the hope that this will be 

O ' measured in future experiments. The predicted results have a rich structure - in particular we find 

a dimensional crossover which can be tuned by varying either the width of the wire or its resistance 
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I. INTRODUCTION 



Recent experiments on disordered superconducting wiresocl have led to renewed interest in the interplay of 
interaction and disorder in low-dimensional systems. These experiments have been made possible by the development 
of fabrication techniques using either stencils or carbon nanotubes, and involve wires of width as low as 50A and 
thickness as low as lOA. An obvious question to ask in such systems is how the transition temperature is affected by 
^ wire width and thickness. Decreasing the thickness of the wire increases the normal state resistance per square which 
is a measure of the amount of disorder in the system. This increase in disorder reduces the transition temperature 
by causing the electrons to move diffusively rather than ballistically, as discussed below. Decreasing the width 



of the wire leads to a crossover from two- to one-dimensional behavior, where increased fluctuation effects further 
• reduce the transition temperature. Ultimately the combined disotder and fluctuation effects completely suppress 
fT^ ■ superconductivity, and a superconductor-insulator transition resultaj. The effect of dimensionality on this quantum 
: phase transition is obviously very interesting from a theoretical point of view. As a first step to understanding this 
transition, we will compare predictions of transition temperature to experimental results. Unfortunately there is very 
little data available for comparison because each data point requires the fabrication of a new wire which is very costly 
in terms of experimental effort. We therefore suggest the application of a perpendicular magnetic field as a way to 
generate more data points for ffifii wire and make predictions for this effect. 
■ The most basic model usedLTu consists of electrons moving in a disordered background interacting via the Coulomb 

repulsion and a featureless phonon-mediated attraction. Simple as this model is, signiflcant difficulties arise even when 
^ . calculating the suppression of transition temperature in two-dimensional films. These difficulties will be compounded 
'■^ ' when looking at quasi one-dimensional wires, which are in a crossover regime, so it is most important that they 
Ch . are fully understood in the simpler two-dimensional case. The basic physics involved is not in dispute, and can be 
O ■ understood from the BCS-MacMillan formula for the transition temperature, Tc, 



T.^1.13..cxp(-^^^^^^1— (1) 



■ where lujj is the Debye frequency (a typical phonon frequency), N{0) is the electronic density of states per spin at 
the Fermi surface, 7 is the featureless attraction, and /i* is the Coulomb pseudopotential, a measure of the Coulomb 
repulsion. The disorder causes the electrons to move diffusively rather than ballistically making them less efficient at 
screening the Coulomb repulsion. This increases /i* and decreases N{0), both effects leading to a decrease in Tc- We 
shall assume that the disorder does not affect phonon properties such as wn-and 7, an assumption that is backed up 
by detailed calculation on at least one model of electron-phonon interactionlJ. Other authorsEJ have questioned this 
assumption, and we intend to address this point in future work. 

The first difficulty which occurs in two dimensions relates to the low-momentum singularity of the disorder- 
screened Coulomb potential, which is given by 
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Vc(q,^m)- Dq^^^ + Dq^ + \u;„,\'' 2N{0) Dq^ ' 

where qxp is the Thomas- Fermi screening wavenumber defined by q^p — 87re^A^(0), D ~ VpT/2 is the diffusion 
constant, vp is the Fermi veiocity,|^d r is the eiastic scattering time. This singuiarity ieads_to a strong suppression 
of the one-particle density of statesE3 which is seen experimentally in tunneling measurementaij. Naively inserting this 
density of states into Eqn. [Q predicts, a large suppression of Tc which is not seen experimentally. This prediction is, 
in fact, incorrect due to a cancellatioio'Ll between the various diagrammatic-Gontributions in perturbation theory. The 
apparently mysterious cancellation turns out to be due to gauge invariancellil and is thus expected to be a very robust 
phenomenon. The end result is that the Coulomb interaction is effectively featureless with strength N{0)Vc = 1/2. 
The fact that we can legitimately use a featureless interaction makes calculations much simpler. Our first task is 
therefore to check that this replacement is still valid in quasi-lD systems. 

The other difficulty is that perturbation theory cannot take us to the superconductor-insulator transition as it 
is only valid for weak disorder, j-and ad hoc attempts to extend its range of validity lead to unphysical reentrance 
problems. Oreg and Finkel'steint^ overcame this problem by developing a resummation technique, which has recently 
been shown to be related to strong-coupling theoryll3. We reexamine their predictions for wires, and find, rather 
surprisingly, that the problem is extremely sensitive to the boundary conditions used. This unexpected sensitivity 
leads us to check all details and approximations used in the resummation technique. Upon satisfactorily completing 
the necessary checks, we compare theory and experiment. Although good agreement is obtained when we focus on 
the suppression of Tc due to increasing disorder, it is less good when we focus on the additional suppression due to 
decreasing wire width. This is largely due to the limited amount of experimental data (only about four data points 
for each of four wire widths) . We propose the application of a perpendicular magnetic field as a means of generating 
more data with a minimum of effort. We extend the resummation technique to include a magnetic field, and make 
predictions for transition temperature as a function of magnetic field for different wire widths and thicknesses. The 
most interesting feature which emerges is a crossover from two- to one-dimensional behavior which can be tuned by 
either decreasing wire width or increasing the resistance per square. We hope that the appropriate experiments will 
be carried out in the near future. 

The remainder of the paper is organised in the following manner. In section (II) we consider first-order perturba- 
tion theory in the absence of a magnetic field. This is the simplest approach possible, and we will use it as a testing 
ground for the validity of the approximations used. Specifically we will show that it is legitimate to use a featureless 
Coulomb interaction and ignore self-energy effects, but that the details of boundary conditions are important. In 
section (III) we make predictions for the transition temperature using the resummation technique, all the necessary 
approximations having been justified in section (II), and then compare theory and experiment. In section (IV) we 
consider the effect of an applied magnetic field on the transitipa temperature. The dimensions of the wires are such 
that we must consider both orbital and Zeeman (spin) effectsc2l. As yet there is no experimental data available for 
comparison to the theory. 



II. FIRST ORDER PERTURBATION THEORY 



We will now calculate the transition temperature, Tc, which we identify as the temperature at which the pair 
propagator, L{q,il = 0), diverges. The latter is defined in Fig. 1(a) in terms of the BCS interaction, 7, and the pair 
polarization function, P(g, 0), by 



L(g,0)-i=7-'-P(g,0). 



(3) 



Let us first consider mean field theory, which corresponds to zcroth order in disorder perturbation theory, for which 
the expansion parameter is l/kp£, where kp is the Fermi wavenumber and £ — vpT is the elastic mean free path. The 
mean field pair polarization function is shown in Fig. 1(b), and leads to the pair propagator 



Lo(g,0)-i =A^(0) 
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The first-order perturbative corrections to the pair polarization function, (5P(g, 0), are shown in Fig. 1(c), and lead 
to a suppression in transition temperature 
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Note that we evaluate the pair propagator at non-zero external momentum, q, since we will eventually consider the 
effect of a magnetic field. 

The five diagrams of Fig. 1(c) can be understood in physical terms: diagrams 1 and 2 are self-energy diagrams 
and give corrections to the electronic density of states, A^(0); diagrams 3 and 4 are pseudopotential diagrams and 
give corrections to the effective Coulomb repulsion, /i*; diagram 5 involves the effect of Coulomb interaction on 
superconducting fluctuations. We will first perform the calculations exactly, and demonstrate the cancellation of the 
low-momentum singularity in the Coulomb interaction which occurs between diagrams 1-4 and 5. We show that this 
is equivalent to considering only diagrams 1-4 with a featureless Coulomb interaction of strength N{0)Vc = 1/2. 
In other words, the only purpose of diagram 5 is to remove the low-momentum singularity, as required by gauge 
invariance. Finally we show that the self-energy diagrams 1 and 2 can be ignored as they give a much smaller 
contribution than diagrams 3 and 4. The contributions from the five diagrams are as follows: 
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aL ^ [Dq^ + 2\ei\f{D{q' + qf + |e,| + |e, +cu,„|) 



i^Vc{q' \Lo„i\)e{-ei{ei + Lo,n)) 
Vc{q' ,\uJrn\)0{ei{ei + LO,n)) 



xp /vmirVrV — V -t- + + +^m| 

3 ^ luj 2^ 2^ aL^ {Dq^ + 2\ei\){Dq^+2\ei+u:^\){Dq'^ + \ei\ + \ei+LU^\y 



5P. = -^-N{Q)tY,tY,^Y. 
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aL ^ [Dq^ + 2\ei\){Dq^ + 2|e, + uJ„,\){D{q' + qf + |e,| + |q + 
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where e/ — {21 -\- \)'itT is a Fermi Matsubara frequency, tOm = 2m7rT is a Bose Matsubara frequency, and a, L are 
the width and length of the wire respectively. The diffusive form of the denominators in Eqn. is valid only for 
momentum and frequency values such that Dq'^, |e/|, ItUml < 1/t, so we introduce the upper cutoff, l/r, in all these 
sums. 

For the rest of this section we will work at zero magnetic field, so we set q — 0. The terms (5Pi_5 have previously 
been summed for 2D filmsu to give the expression 
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The 1/q'^ singularity in the screened Coulomb potential Vc{q' , ojm) is cancelled since there is an overall factor of Dq'^ 
in the above equation. Note that there is, in fact, no divergence at Dq'^ = ujm, but this can only be seen by Taylor 
expanding to second order. The need to perform this expansion makes the calculation involved, especially at the level 
of numerics. To allow simplification, we define y = Dq'^ /2'kT and m = uim/2iTT, which leads to the result 
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where we have spht off the most divergent term at large momentum and frequency, and M — 1/2'hTt. To proceed 
further we must consider the q'-sum in various dimensions. In two dimensions the sum becomes an integral, and it is 
most sensible to perform the cut-off in a circularly symmetric manner 
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In quasi-lD wires, we must replace the integral in the transverse direction by a sum, 
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where q'^ ~ q'^ + q'^ ■ At this point we have to decide which boundary conditions to use in the transverse direction. To 
illustrate the strong effect this has we will compare periodic and zero-gradient boundary conditions, the details of which 
are listed in Table 1. The zero-gradient Jafiundary conditions are the physically correct ones for a superconductor- 
insulator boundary as discussed in Ref. To perform the sum over q' note that we can write y — Ul + Ua where 
yL — Dq'l/2'KT and ya — Dq'^ /2t:T, and if we define z — ^JyZ, the g^-integral becomes 
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(11) 



where Lt = \J2'kD jT is the thermal length. Note that y = + Ap{a/LTY for zero-gradient boundary conditions 



and y 



■ p^a/Lx) for periodic boundary conditions. The final result of perturbation theory is then given by 



Eqn. (H) with the prefactor and sums becoming 
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In the above, t is a dimensionless measure of disorder given by 
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\Q2kQ,. Note that t only differs from the 



where Ru is the normal state resistance per square, and Rq = 2TTh/e^ 
previously discussed parameter l/kp£ by a factor 27r. 

Having shown that the low-momentum singularity in the screened Coulomb potential is cancelled by diagram 5, 
we recalculate the suppression of transition temperature using a featureless Coulomb potential of magnitude A^(0) Vc — 
1/2, and including only diagrams 1-4. We have previously showr£3 that the above equations simplify to give 



SPl2 



SP., 



-^E 

^2TE 



2ttT 



) 



2'kT 



E (£)q/2 

q' 



(14) 



For convenience we define dimensionless sums Ii (m) and I2 (m) as 
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which are evaluated in two dimensions and quasi-one dimension in Appendix A. When we combine Eqns. (|^) 
( |l^ ) we obtain the formula for transition temperature suppression assuming a featureless Coulomb potential, 
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If we ignore the self-energy diagrams (^^12), the him) term in the above equation vanishes to leave 
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Having derived expressions for the suppression of transition temperature using various approximations and bound- 
ary conditions, we perform a detailed numerical analysis, the results of which can be summarised as follows. The 
transition temperature is very sensitive to the boundary conditions used - periodic boundary conditions do not yield 
a dependence of transition temperature on wire width, whilst zero-gradient boundary conditions show a systematic 
suppression as wire width is decreased. The predictions are also sensitive to the value of the thermal length, Lt, since 
the width, a, enters only in the dimensionless combination a/ Lt- Experimentally Lt may be known only roughly and 
it may be necessary to fit its value. On the other hand, the predictions are not very sensitive to either the form of the 
Coulomb interaction used (Eqns. (||) and (p^), or the inclusion of self-energy corrections (Eqns. ( p^ ) and (p7|)). The 
slight difference between these approaches vanishes when we fit the initial slope of the TdRa) curve to experiment, 
which sets the value of the BCS attraction, 7. This demonstrates that the simplest approach - featureless Coulomb 
potential and no self-energy diagrams - is valid. This is important since these simplifications are necessary for the 
resummation technique discussed in the next section. It is necessary to use this technique to cure the unphysical 
re-entrance problems which beset perturbation theory at low temperature. 



III. THE RESUMMATION TECHNIQUE 

In this section we will review the resummation techniqueEl developed by Oreg and Finkel'stein (OF) as applied 
to the problem of suppression in quasi-lD wires. This method involves calculating the pair-amplitude, r(e„,ei), 
which is a matrix with Fermi Matsubara frequencies labelling its rows and columns. The transition temperature is 
identified as the temperature at which this matrix becomes singular. OF obtained an equation for F by including ouly 
pseudopotential diagrams, and solved it using matrix perturbation theory with zero-gradient boundary conditionsEZI. 
We have performed calculations to show that it is legitimate to include only these diagrams, and that it is essential 
to use zero-gradient boundary conditions, as described in section II for perturbation theory. In a previous papeillj we 
have shown that the equation for F can be solved exactly at no extra numerical cost, so the matrix perturbation theory 
is unnecessary, and indeed gives incorrect results at large disorder. Therefore we shall include only pseudopotential 
diagrams and solve the equation for F using an exact diagonalization method with zero- gradient boundary conditions. 

Before we proceed further, let us briefly describe the details of the technique. The ladder summation involved 
is shown in Fig. 2. The block iA, an effective Coulomb pseudopotential, contains diagrams 3 and 4 of Fig. 1(c); 
the block S, the self-energy of the Cooperon impurity ladder, C, contains diagrams 1 and 2 of Fig. 1(c). We have 
included the self-energy diagrams for completeness, although their effects can be ignored in calculations, as discussed 
above. The equation for F is seen to be 

M ^ 

F(e„, ei) = -I7I + <A(e„, ei) ~ 2ttT ^ [-I7I + tA(e„, e,„)] _F(e™, e,) (18) 

m— 

where the contributions to tA are 

2^iV(0)aL^ [W^ + |6„| + |e™|P ^ '"'"^ 

tA4ien,e,n)^ A^^n^ ~^E rn I 1'^ I I I — rr^(^»'^™)- (19) 
27riV(0) aL ^ [Dq'- + |e„| + |e,„|J 

We have kept non-zero external momentum q so that we can include magnetic field effects later. For now we consider 
the case of zero magnetic field, so we set q = 0, and the two contributions to tA are identical. It follows that A(e„, e™) 
does not depend upon the relative sign of e„ and Cm and has the form 
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A(e„,e„0 = + 1), (20) 
where Ii is defined in Eqn. (^5|). Eqn. (|l^ ) can then be written in matrix form as 

f = -|7|i + tA-[-|7|i + <A]f:-if, (21) 

where f = r(e„, €,„), = 1, A„„i = A(e„, €„) and e„m = (n + l/2)5nm- This has the solution 

f = i'/'ii - |7|n)-ie-i/2(-|7|i + tA), (22) 

where 

n = e-i/2[i - |7|-4A]e^i/2^ (23) 

and /„,„ — 6nm is the identity matrix. It foUows that when the matrix 11 has an eigenvalue equal to I/I7I, the pair 
amplitude diverges, and we have found Tc- In other words, we need to solve the eigenvalue equation 

|7r^/-n(r,)] iv-) = 0. (24) 

We solve Eqn. ( p^ ) by diagonalizing the matrix [|7|^^/— n(r)] numerically, and decreasing the temperature 
until its lowest eigenvalue equals zero. Since we only need the lowest eigenvalue, and the matrix has such a simple 
structure, the Lanczos method is very efficient here. The only minor complication is that the matrix n(r) depends 
upon T both through the dependence of A upon T and through its rank M — 1/2ttTt. We therefore start at the 
mean- field value of M, which we will call Mq, and decrease the temperature by increasing M successively by one. We 
diagonalize the matrix [|7|^^/ — II] for each value of M until the lowest eigenvalue changes sign. At this point we have 
found Tc for the given problem, and Tc/Tco — Mq/M. This is obviously the method of choice, since no approximation 
is involved. 

We can now compare thftory and experiment: rihere are three sets of esperimental data that we know about - 
those of Graybeal et al (1987)y, Sharifi et al (1993)0, and Xiong et al (1997)0. The data of Graybeal et al consists of 
Tc vs a for wires of constant thickness (and hence constant resistance per square) . Their data shows that the relative 
suppression of Tc varies as 1/a^ with a proportionality constant independent of R^. We have performed calculations 
which predict a dependence roughly of the form 1/a, but with a proportionality constant dependent on i?^. We will 
therefore not consider this data further. The data of Sharifi ct al and Xiong et al are both of the form of Tc{Ra, a) vs 
Ra for wires of different width, but only the latter shows a consistent trend with decreasing width, so it is only this 
data that we shall attempt to fit. The fact that this is the most recent data is reassuring, since it suggests that more 
recent fabrication techniques have produced better quality wires. 

In Fig. 3 we compare the data of Xiong et al with our theoretical predictions. In doing this we encounter 
two difficulties. The first is that experiment and theory do not agree for some of the 2D films used. We assume 
that this is due to some additional physical processes that will affect samples in the same manner independent of 
width. We therefore multiply our theoretical predictions by Tc(2D,expt)/rc(2D, theory) to factor out this effect. The 
second difficulty is that wires of different widths do not appear to have the same value of Tc in the clean limit, 
Tc{0, a). We therefore assume a different Tc(0, a) for each wire. Upon making these two assumptions, we see that the 
agreement between theory and experiment is good. However plotting the data in this form emphasises the suppression 
of transition temperature due to disorder, whereas we are interested in the additional suppression due to the decrease 
in width. For this reason we need to replot the data to emphasise this latter point. 

In Fig. 4 we replot Fig. 3 in terms of the additional relative suppression, 0(i?n, a), which is defined by 

n/D ^ TciB„,a)/Tc{0,a) 

^^'^'''^ = Tc{R.,2D)/Tc{0,2Dy ^^'^ 

This quantity focuses on the basic physical phenomenon we are interested in - the additional suppression of transition 
temperature caused by the reduced dimensionality in the quasi-lD wires. It has the additional benefit of factoring out 
the two difficulties discussed in the previous paragraph. When plotted in this form, the agreement between theory 
and experiment no longer appears to be as good. This is because our theory, and indeed any theory, predicts a smooth 
dimensional crossover whereas the experimental data shows an abrupt crossover between the 2D film and the lOOOA 
wire, which lie on one curve, and the 580^1, 350A, 250A wires, which lie on a different curve. The fact that the data 
for the three narrowest wires fall on top of each other might be due to them having the same effective dimensionless 
width a/Lx- If we allowed ourselves the additional freedom of fitting this parameter separately for each wire, we 
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would obtain good agreement between theory and experiment, but we do not believe this is justified. We have used 
the value of the thermal length quoted in Sharifi et al (viz Lt = 400A at T = AK, which corresponds to Lt = 300A 
at Tco = &.%K). 

The basic problem is clearly that there is not enough experimental data available - there are only 19 points in 
total, and each of these corresponds to a different wire. To carry out a more systematic comparison between theory 
and experiment, one would need wires of many different resistances (thicknesses) for each value of width. This would, 
however, be very costly in terms of experimental effort. A less costly method for generating more data would be to 
apply a magnetic field to a given wire. We therefore extend our theory to include the effect of such a field in the next 
section. 



IV. EFFECT OF AN APPLIED MAGNETIC FIELD 



We will now consider the effect of an applied perpendicular magnetic field. There are two effects to consider here: 
the orbital effect, caused by the magnetic field scrambling the relative phase of the electrons in a Cooper pair, and 
the Zeeman effect, caused by the magnetic field lifting the degeneracy of up and down spins. Let us first consider the 
relative importance of the two effects. The orbital effect will destroy superconductivity when the magnetic dephasing 
rate, 1/th, is of the order of ksTc, which gives H w AT for a Pb wire of width 250 A. The Zeeman splitting will 
destroy superconductivity when its magnitude (J-bH is of the order of ksT^, which corresponds to ~ lOT for Pb. 
Since these fields have the same order of magnitude, we will need to consider both effects. 

We first examine the orbital effect. The magnetic field causes the eigenvalues of the Cooperon operator to change 
from the zero-field values 
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(26) 



(27) 



where we have chosen the gauge A ~ (—Hy, 0,0). The eigenfunctions 'ip{x,y) must satisfy zero-gradient boundary 
conditions on the sides of the wire, y — ±a/2, and periodic boundary conditions at the ends of the wire, x = 0,L. 
The solution has the separable form 



i^{x,y) = 0,„,p(?/)e*'?^, 



2?™ 
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(28) 



where 



,p 



satisfies the equation 



D 



'^qn,p 



dy 



^(±«/2) = 0, 



(29) 



and p = 0, 1, 2 . . . labels the eigenvalues in ascending order. 

To solve this eigenvalue problem numerically, we convert it into a matrix diagonalization problem by working in 
a truncated function space spanned by eigenfunctions of the zero-field problem. The lowest eigenvaluey,Ao,Oj equals 
2/th, where 1/th is the magnetic dephasing rate which appears in the Abrikosov-Gor'kov equationEJ for Tc{H). 
There are two asymptotic limits where 1/th is easily determined: the 2D limit, where the problem becomes that of 
a simple harmonic oscillator, with 1/th = DeH; and the ID limit, where the eigenfunction is simply a constant, and 
hence (using perturbation theory) 1/th — De^ H^a? /Q. The 2D to ID crossover occurs when the width, a, is of the 
order of the magnetic length, Lh — {^/eHy^^, which occurs at the field Hq = h/ea?. At larger magnetic fields the 
cyclotron orbit of the electron lies entirely within the wire, and hence we are in the 2D limit; at smaller magnetic 
fields the electron hits the sides of the wire before completing a cyclotron orbit, and hence we are in the ID limit. To 
summarise, 
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These eigenvalues then enter into our equations for the pair propagator in two ways: (i) the term Dq^ in the 
q-dependent pair propagator is replaced by the lowest eigenvalue Ao,o; (ii) sums over functions of Dq'^ in diagrams 
which involve Cooperonic ladders (e.g. diagrams 2 and 4 in Fig. 1(c)) are replaced by sums over Xq^^p. Obviously 
sums over functions of Dq'^ in diagrams which involve diffusion ladders (e.g. diagrams 1 and 3 in Fig. 1(c)) remain 
unchanged. It is numerically very expensive to perform the sums over all the eigenvalues, so we would like to see how 
important it is to use these eigenvalues in any sums. To answer this question we again use perturbation theory as our 
testing ground. We perform the Cooperonic sum in two different ways - both correctly, and by replacing it by the 
corresponding diffusion sum. In both cases we still make the replacement of Dq^ by Ao,o = '^/th- The perturbative 
equations are then 
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where a = 1/2'itTth- The diffusive sums /i(to) and l2{'nfi) are defined as in Eqn. ( |l5| ) except that here we perform 
the momentum sums in the longitudinal direction exactly to give 
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and the Cooperonic sum, Ci , which is now distinct from Ii , is given by 
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(33) 



We found essentially no difference between using the correct Ci , or replacing it by Ii , over a wide range of magnetic 
fields. It follows that the only replacement we need to make is of Dq^ by Ao,0 7 which greatly reduces the numerical 
cost. 

Let us next consider the effect of the Zeeman splitting. This leads to changes in both the mean field and 
perturbative terms, and we will show that only the former effect is important. To do this we will evaluate the 
mean-field pair polarization function and the first-order perturbation correction, 5P, both with and without Zeeman 
splitting. We find that the latter is essentially unaffected. To treat the Zeeman effect, we need formulas for the 
Cooperon impurity ladder and Hikami boxesEj, where we note that the boxes in diagrams 1 and 3 of Fig. 1(c) will 
give different contributions. The non-zero elements of the spin-dependent Cooperon ladder, Cap-ySi ar6 
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27riV(0)T2 Dg2 



(34) 



where A = /i^iJ is the Zeeman splitting. The appropriate spin combination which occurs in an s-wave pair propagator 
is C|||| + C||.tt = 25fteCf-fjj. This leads to the mean field formula for Tc, 
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5Re 
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(35) 



where n = a — i(3 with /3 = /S./2'kTc = ^bH /"^-t^Tc. To calculate the first-order pair polarization terms ^Pi_4, we note 
that: (a) the diffusion ladders on Coulomb vertices are unchanged; (b) Cooperon ladders pick up an extra =F2iA as 
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described above; (c) the Hikami box in 6Pi picks up an extra T2iA, whilst the box in SPs is unaffected. This leads 
to the formulas 
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where I\{rn) and I^ijn) are the same as before, whilst 
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(37) 



Calculating the 5Pi-4 as above, or with /3 set to zero, we find essentially no difference over a wide range of magnetic 
fields until wc get to widths of order 10 A. It follows that we only need consider the Zeeman splitting at the mean 
field level. Even this is a small effect, only noticeable for wires narrower than about 150A. 

Let us finally perform calculations of using the resummation technique. From the discussions above, we now 
know that to include magnetic field effects we need only replace Dcf' in Coopcrons by ^/th, which must be calculated 
numerically, and include the Zeeman splitting at the mean field level. The Cooperon denominator Dtf + 216^1 then 
becomes 2\em\ + 2/Tij — 2iA, and hence Cnm = {n + 1/2 + n)5nrm which is now complex. Because the two spin 
configurations which occur in an s-wave pair propagator are complex conjugates of each other, we end up taking the 
real part of the complex matrix equation. The equation for A„m becomes 



_ r Ji (n + m + 1 + 2a) 

~ 7i(n + m+l + 2a)+272(n + m+l + 2a) 



enCm > 
e„em < 



(38) 



where we have not included the Zeeman effect as wc have shown this to be much smaller in perturbation theory. 

We are now in a position to calculate the transition temperature, Tc, as a function of resistance per square, 
Rn, wire width, a, and applied magnetic field, H. This obviously leads to a wide choice of ways to present results. 
Experimentally one would probably vary magnetic field on a sample of fixed width and thickness, and measure Tc{H). 
In Fig. 5 wc therefore plot Tc{H) for wires of different resistances per square, Ra = OO, lOOOfl, 2000fi, and widths 
a = 500A, 600A, 750A, 2D film. There is far too much data on this plot, but some general trends are apparent. 
Three physical processes are occuring here: increasing increases the effective Coulomb repulsion which strongly 
suppresses T^, decreasing a means that a larger magnetic field is required to fit a flux quantum inside the wire, and 
hence Hc2 increases at fixed Ra] decreasing a also causes a crossover from two to one-dimensional behavior, and so 
increases fluctuation effects, leading to a small reduction of Tc at fixed Rn- The last two effects mean that the Tc{H) 
curves for wires with the same Ra and different a cross each other at non-zero field. 

Another effect of the dimensional crossover is to change the shape of the Tc{H) curve. This occurs because the 
magnetic pair-breaking rate, 1/th, is linear in magnetic field in the two-dimensional limit, but quadratic in magnetic 
field in the one- dimensional limit. This leads to two limiting forms for the dimensionless Tc{H) curve. 
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(39) 



where p = 1 in the 2D limit and p = 2 in the ID limit. We can move from the 2D limit to the ID limit either by 
decreasing a at fixed as shown in Fig. 6, or by increasing R^ at fixed a, as shown in Fig. 7. Considering first 
the width dependence, we see in Fig. 6 that the dimensional crossover occurs when a ~ Lt- For the wider wires 
the Tc{H) curve joins the two-dimensional curve at the value of the magnetic field such that the magnetic length, 
Lh, is of order the width, a; for narrower wires this value of H exceeds the two-dimensional value of Hc2, and the 
curves never meet. For ultranarrow wires, a < 250A, the Zeeman effect becomes important. This causes the Tc{H) 
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curves to change shape again, moving slowly back in the direction of the 2D limit. Now let us consider the resistance 
dependence of the Tc{H) curve. For narrow wires, a « 250A, all Tc{H) curves fall on the ID curve independent of 
i?n. For wide wires, a ~ lOOOOA, all Tc{H) curves fall on the 2D curve independent of Ra- For wires of intermediate 
width, the Tc{H) curves can be changed from the ID curve to the 2D curve by increasing Ra, as shown in Fig. 7. 
Finally for ultranarrow wires, a « 50 A, the Tc{H) curves can again be tuned by varying R^, in this case the tuning 
being between one-dimensional orbit and Zeeman curves. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have considered the effect of localization and interaction on transition temperature in ultranarrow 
superconducting wires. Our main conclusion is that, taking into account the relatively small amount of experimental 
data, theory and experiment are consistent with each other. The agreement is not perfect: the experimental data 
seems to show a very sudden crossover from 2D to ID behavior when the wire width becomes of the order of the 
thermal length, whilst theory predicts a smooth crossover. Indeed, when plotted as the additional relative suppression 
of transition temperature due to finite width, the experimental data falls on two curves - one for the 2D film_and 
lOOOA wire; the other for the 580A, 350A and 250A wires. These conclusions are based on a single experimenta, as 
this is the only experiment which shows any consistent trend with wire width. More experimental data is obviously 
needed to investigate this further. Given the difhculty in fabricating these wires, we suggest that the magnetic field 
dependence of transition temperature be measured in future experiments. We have therefore made predictions for 
magnetic field dependence of transition temperature in wires of varying width and resistance in the hope that such 
measurements will be made in the future. Our main prediction is of a change in the shape of the reduced Tc{H) 
curves due to the dimensional crossover. This may be tuned either by varying wire width-or resistance. 

Our work is based on the resummation technique developed by Oreg and Finkel'steinEj, which is a very powerful 
method of going beyond perturbation theory, and hence giving believable results at strong disorder. The suppression 
of the transition temperature is very sensitive to boundary conditions, due to the different relative weights of the 
zero momentum mode. We have used the physically correct zero-gradient boundary conditions in our work. The 
unexpected strong dependence of results on boundary conditions led us to worry about all details of the calculation. 
We therefore also investigated the effects of using a featureless Coulomb interaction instead of the exact one, ignoring 
self-energy contributions, and solving the OF matrix equation exactly instead of using an approximate method. We 
found that it is legitimate to use a featureless Coulomb interaction and ignore self-energy contributions, but that our 
matrix equations should be solved by exact diagonalization (which in any case is no more difficult). 

In our calculations of the magnetic field dependence of transition temperature, we find that the orbital effect is 
the most important, and can be treated by replacing the Dq'^ terms in the momentum dependent pair propagator, 
L{q,0), by twice the magnetic dephasing rate, 2/th- The latter is the lowest eigenvalue of the minimally coupled 
Cooperon operator, and must be found numerically. The Zeeman splitting is found to have a non-negligible effect 
only at the mean field level, and even this effect is small. The perturbative term only becomes important in very 
narrow wires (of order lOA), which cannot currently be fabricated. 
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APPENDIX A: DETAILS OF SOME MOMENTUM SUMS 



In this appendix we derive expressions for the quantitities Ii (m) and I2 (m) defined in Eqn. ( [iq ) as 
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aL ^ iDq'^ + \u^m\y ^ ' 
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In two dimensions the sum becomes an integral, and cutting off in a circularly symmetric manner, 
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we obtain the results 
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In quasi-lD wires we replace the integral in the transverse direction by a sum, so that /i(m) becomes 
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2TrTm 
D 



: tan 
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(A4) 



and him) can be found by differentiating the above with respect to m and dividing by —2'kT. Substituting in the 
two types of boundary condition listed in Table 1 yields 
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and 
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where A = a/Z/T is a dimensionless width, the ratio of the wire width to the thermal length. 
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Type 


Definition 


Momenta, qa 


Details of sum over p 


Periodic 


ip{x + a) = tp{x) 


2'np 
a 


PO 

E 

P=-PO 


Zero- Gradient 


tp'{x = ±a/2) = 


■np 
a 


E 

p=0 



TABLE I. Details of the boundary conditions discussed in the text. The upper cut-off in momentum is given by Dq^ = 1/r, 
so that po = a/2TT^ Dt = ^fM{a/ Lt) where Lt = \/2'kD/T is the thermal length and M = 1/2'kTt. 
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(la) (lb) (Ic) (2) 




(c) 

FIG. 1. The diagrams needed to calculate the first-order perturbative correction to the pair-propagator L(q, 0). (a) Definition of 
pair propagator, L{q, 0), in terms of the pair polarization function, P{q, 0), and the 4-point BCS interaction — 17|. (b) Zeroth- 

ordcr (mean-field) polarization function, Pq (g, 0) , in a dirty superconductor. The solid linos arc electron Green functions; the 
dashed lines are impurity lines, (c) First-order perturbative correction to pair polarization function, 6P{q, 0). The wiggly line 
is the screened Coulomb interaction. 
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FIG. 2. Diagrammatic equation for the pair amplitude matrix r(e„,ei). Block 7 is the BCS interaction. Block tA is the 

correction to the cfTective interaction caused by the interplay of Coulomb interaction and disorder. Block S is the correction 
to the bare Cooperon ladder caused by Coulomb interaction and disorder, which yields a negligible contribution and is thus 
ignored in our numerical calculations. 
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I ~ \ ~ \ ~ 1 

1000 2000 3000 

Resistance R (Q.) 

FIG. 3. Comparison of theory (solid lines) and experiment (symbols) for transition temperature, Tc, as a function of resistance 
per square, Ra, for various wire widths, a — 250A, 350A, 580A, lOOOA. The experimental data is that of Xiong et afl. The 
theory used is the one we have argued to be correct - resummation technique solved by exact diagonalization using zero-gradient 
boundary conditions. The theory plots are, from top to bottom, for lOOOA, 580A 350A and 250A wires. Rescaling has been 
applied to take into account the fact that theory and experiment do not agree with each other in two dimensions, and that the 
transition temperature at zero resistance appears to depend upon midth. The thermal length, Lt, has been taken to be 300A 
at Tco = 6.9K throughout - this is the value quoted by Sharifi et al3. Plotting the data in this form focuses on the suppression 
of transition temperature due to disorder, and we see that the agreement between theory and experiment is good, particularly 
for the lOOOl and 2501 wires. 
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1.2 




FIG. 4. Comparison of theory (solid lines) and experiment (symbols) for the additional relative suppression of transition 
temperature, &{Ra,a), as a function of resistance per squaxe, for various wire widths, a = 250A, 350A, 580A, lOOOA. 
Plotting the data in this form focuses on the suppression of transition temperature due to width, and we see that the agreement 
between theory and experiment is not particularly good. This may have been anticipated from the experimental data which 
shows an abrupt crossover between the 2D film and the lOOOA wire, which lie on one curve, and the 580A, 350J1., 250A wires, 
which lie on a different curve. Our theory, and indeed any theory, predicts a smooth crossover. 
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FIG. 5. Predicted transition temperature, Tc, as a function of applied perpendicular magnetic field, H, for wires of various 
widths, a = 500A, 600A, 750A, 2D film, and resistances per square, = Ofl, lOOOfi, 2000ri. Three physical effects can be 
seen in this figure: increasing reduces Tc dramatically due to the increased Coulomb repulsion; decreasing a increases the 
upper critical field, i?c2i as a larger field is needed to fit a flux quantum into the wire; decreasing a at fixed Ra causes a small 
additional suppression of at zero field due to increased fluctuation effects caused by the crossover from two to onc-dimcnsional 
behavior. The last two effects mean that the Tc{H) curves for wires with the same and different a cross each other at 
non-zero field. Finally the shape of the Tc{H) curve changes as width a is altered due to the dimensional crossover. There is 
obviously a large amount of data in this figure, and to analyse the individual trends we will need to re-plot some of this data 
in a different form. 
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FIG. 6. Predicted dependence of the shape of the Tc{H) curve on wire width, a, at a fixed resistance per square, Rn — lOOOO. 
We plot using dimensionless variables Tc{H)/Tc{H = 0) vs H/Hc2 in order to compare the curve shape. As width is reduced 
we see a crossover from the two-dimensional limit, where the magnetic pair-breaking rate 1/th oc H, to the one-dimensional 
limit where 1/tjj (X H^. The crossover occurs when the wire width is of the order of the thermal length, L^- For the wider 
wires the Tc{H) curve joins the two-dimensional curve at the value of the magnetic field such that the magnetic length, Lfj, 
is of order the width, a; for narrower wires this value of H exceeds the two-dimensional value of Hc2, and the curves never 
meet. For ultranarrow wires, a < 250A, the Zeeman effect becomes important. This causes the Tc{H) curves to change shape 
again, moving slowly back in the direction of the 2D limit. 
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FIG. 7. Predicted dependence of the shape of the Tc{H) curve on resistance per square, Bu, at a fixed wire width, a = 1000 A. 
We see a crossover from two-dimensional to one-dimensioanl behavior as is increased. This disorder-tuned crossover occurs 
only for wires of intermediate width; a similar plot for a = lOOOOA shows all data falling on the two-dimensional curve 
independent of i?n, whilst that for a = 250A shows all data falling on the one-dimensional curve. 
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